Evidence for a long-range component in the pion emission source in Au+Au collisions 

at v /s^7=200 GeV 
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Emission source functions are extracted from correlation functions constructed from charged pions 
produced at mid-rapidity in Au+Au collisions at % /Sjvjv=200 GeV. The source parameters extracted 
from these functions at low kr, give first indications of a long tail for the pion emission source. The 
source extension cannot be explained solely by simple kinematic considerations. The possible role 
of a halo of secondary pions from resonance emissions is explored. 

PACS numbers: 25.75.Ld,25.75.Dw 



Collisions between ultra-relativistic heavy ions can The decay dynamics of this matter are strongly influ- 
lead to extremely high energy-density nuclear matter . enced by the nuclear Equation of State (EOS) and possi- 
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bly by a de-confined phase ■ An emitting system which 
undergoes a strong first order phase transition is expected 
to show a much larger space-time extent than would be 
expected if the system remained in the hadronic phase 
throughout the collision process 0. Indeed, several hy- 
drodynamical calculations show such an increase for par- 
ticle emission sources Q.Q • providing hadronization does 
not occur via a supercooled state |5j. It has also been 
suggested that the shape of the emission source function 
can provide signals for a second order phase transition 
and whether or not particle emission occurs near to the 
critical end point in the QCD phase diagram 

Interferometry studies provide important information 
on the emission source function for particles produced 
in nuclear reactions ranging from elementary collisions 
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and (p)pp) to those involving very heavy ions 
. Recent studies include the use of two- and three-pion 
interferometry correlations spanning the beam energies of 
the Alternating Gradient Synchrotron (AGS), the Super 
Proton Synchrotron (SPS) and the Relativistic Heavy Ion 
Collider (RHIC) {y/hfN~ 2 - 200 GeV) [HI 03, El G3 
A common theme for these papers is the extraction 
of the widths of emission source functions which are 
assumed to be Gaussian. For such extractions, the 
Coulomb effects on the correlation function are usually 
assumed to be separable 01 as weu - Such an approach 
was followed in an earlier analysis in which we used the 
Bowler-Sinyukov 3D HBT analysis method [in Bertsch- 
Pratt coordinates] to probe for long-range emissions from 
a possible long-lived source 0. The RMS-widths so 
obtained for each dimension of the source i?i on g, -Rside 
and -Rout, gave no evidence for such emissions, suggest- 
ing that the sound speed is not zero during an extended 
hadronization period. 

In this letter we exploit the mo del- indep endent imag- 
ing technique of Brown and Danielewicz [161 Il7j ] to make 
a more detailed study of both the shape and the space- 
time characteristics of the pion emission source func- 
tion. The method uses a numerical calculation of the 
two particle wave function (including final state interac- 
tions (FSI)) to produce an inversion matrix that operates 
on the correlation function to produce the corresponding 
source function. The technique has been used to address 
only a few data sets at relativistic beam energies. 

The measurements were made with the PHENIX de- 
tector 20) at RHIC. The charged pions relevant to 
this analysis were detected in the two central arms 
(M < 0.35). Track reconstruction was accomplished 
via pattern recognition using the drift chamber (DC) 
followed by two layers of multi-wire proportional cham- 
bers with pad readout located at radii of 2, 2.5, and 
5 m |20j. Particle momenta were measured with the 
resolution <5p/p = 0.7% ®1.0%p (GeV/c). Very good 
pion identification (PID) was achieved with a 2a cut 
about the pion peak in the squared-mass distribution 
for pr 2.0 GeV/c and pt 1 GeV/c in the Time 



of Flight (TOF) and the Electromagnetic Calorimeter 
(EMC) respectively. The event centrality was determined 
using the PHENIX beam-beam counter and the zero de- 
gree calorimeter |2l|. Approximately 22 million Au+Au 
events were analyzed to generate and study correlation 
functions for several centrality and pr selections. 

Two-pion interferometry correlations were obtained 
via the correlation function C(q) — N cor (q) / 'N m i X {q), 
where the numerator is the relative momentum distri- 
bution of particle pairs from the same event (foreground 
pairs) and the denominator is the relative momentum 
distribution of particle pairs obtained from mixed events 
(background pairs). Here, q = hy~ Jjpi ~ P2) 2 is half 
of the relative momentum between the two particles in 
the pair cm. frame (PCMS). p\ and P2 are the momen- 
tum 4- vectors of each particle in the pair and C(q) = 1 
for large g-values where final-state interactions are neg- 
ligible. Track-pair cuts similar to those of Ref. |l2( 
were applied to foreground and background pairs respec- 
tively. That is, pairs within 0-5 cm in the beam direction 
(AZdc) an d 0.02 radians in azimuthal angle (A(f>oc) m 
the DC were eliminated from the pair sample to remove 
ghost tracks, and pairs within 0.0 < A4>dc < 0.1 radians 
for AZpc > 5 cm were removed to avoid an inefficient 
region. The latter set of cuts were supplemented with the 
removal of pairs having a separation AR < 14 cm and 
< 16 cm in the TOF and EMC respectively. Systematic 
variations of all of these cuts were explored to determine 
systematic error estimates; very little influence on the ex- 
tracted correlation functions was found. Careful studies 
of the influence of momentum resolution on the corre- 
lation function were also made following the simulation 
technique outlined in Ref. |l2^ . The maximum effect was 
found to be ~ 0.4% change in the correlation function at 
low q values, so it was neglected. 

The filled circles in Fig.^a) show the one-dimensional 
(ID) correlation function C(q), for a centrality of - 
20% and for 0.20 < k T = \(p T ,i + Pr, 2 ) < 0.36 GeV/c. 
The characteristic enhancement for q < 25 MeV/c re- 
flects a combination of Bose-Einstein statistics and the 
FSI between pion pairs. The correlation function is not 
Coulomb corrected because the FSI (including Coulomb 
but no strong interaction) is included in the imaging and 
fitting procedure as described below. 

The ID correlation function and source function S{r) 
(the probability of emitting a pair of particles at a 
separation r in the PCMS frame), are related via the 
ID Koonin-Pratt equation |22^ : 



C(q) - 1 =4?r / drr 2 K (q,r)S(r). 



(1) 



The angle-averaged kernel K (q,r) encodes the FSI and 
is given in terms of the final state wave function $ q (r), 
as K (q 7 r) — 5 / d(cos(# qi r)Hj$q(r)| 2 — 1), where 6^ r is 
the angle between q and r [TjJ. The procedure for the 
inversion of Eq.nto obtain S(r) is also given in Ref. [l7| . 
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FIG. 1: (color online) (upper, a) (filled circles) Correlation 
function, C(q) for tt + tt + and iv~ir~ pairs; (open squares) re- 
stored correlation from imaging technique; (dotted line) direct 
correlation fitting; (open triangles) ID angle- averaged correla- 
tion of 3D correlation function, (lower) ID source function (b) 
S(r) and (c) 47rr 2 S(r): (open squares) imaging; (filled circles) 
spheroid fit to correlation function; (open triangles) angle- 
averaging of 3D-Gaussian source function. Systematic errors 
are less than size of data points. 



The open squares in Fig. ^b) show the source func- 
tion obtained from the correlation function presented in 
Fig. ^a). As a cross check of the imaging procedure, 
a restored correlation function was generated via Eq. ^ 
with the extracted source function as input. The open 
squares and filled circles in Fig. ^a) indicate excellent 
consistency between the measured and restored correla- 
tion functions. The ID-source function (cf. Fig. ^b)) 
points to a Gaussian-like pattern at small r and a pre- 
viously unresolved "tail" at large r. The robustness of 
this tail was established via an extensive array of tests 
including its dependence on pair and PID cuts, and on 
momentum resolution; no variation outside of the stated 
error bars was found. 

This new observation of a tail is made more trans- 
parent via a comparison with the source function con- 



structed from the parameters (i?i ng, -Rsid&L-Rout and A), 
obtained in an earlier 3D HBT analysis [l5j. The proce- 
dures outlined in Ref.0] were employed to construct this 
source function (see Fig.^l. The measured and 3D angle- 
averaged correlation functions differ for q < 15 MeV/c as 
do the respective source functions for r > 17 fm. The im- 
aged source function exhibits a more prominent tail than 
the angle averaged 3D HBT source function. This dif- 
ference could stem from the Gaussian shape assumption 
employed in the 3D HBT analysis. The 3D Gaussian fit- 
ting procedure by construction is sensitive only to the 
main component of S(r), and thus would not be capable 
of resolving fine structure at small-q/large-r. Given the 
fact that the volume clement increases quadratically with 
pair separation, this difference is considerable as shown 
in Fig^Jc), where the radial probabilities (47rr 2 S(r)) are 
compared. The open triangles in Fig. ^a) clearly indi- 
cate that the differences in the source functions reflect 
an important disparity in the corresponding correlation 
functions for q < 10 - 15 MeV. 

Parameters of the source function for different as- 
sumed shapes were extracted via direct fits to the cor- 
relation function. Filled circles in Figs. ^b,c) show the 
source function obtained from such a fit for a Spheroidal 
shape [13 , 
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where R c s =1/ (l/R T — l/-Ro), Rt is the radius of 
the Spheroid in two perpendicular spatial dimensions and 
Rq = ax Rt is the radius in the third spatial dimension; a 
is a scale factor. The long axis of the Spheroid is assumed 
to be oriented in the out direction of the Bertsch-Pratt 
coordinate system. The fraction of pion pairs which con- 
tribute to the source A, is given by the integral of the 
normalized source function over the full range of r. 

The procedure for making a direct fit to the correlation 
function involves the determination of a set of values for 
the Spheroid parameters of Eq. [21 which reproduce the 
observed correlation function when the resulting source 
function is inserted into Eq.^ The minimization package 
MINUIT was used to minimize the \ 2 between the ob- 
served and calculated correlation function. The x 2 / n( Af 
value so obtained was ~ 1. The dotted curve in Fig.QJa) 
shows the fit to the data. The resulting source function 
shown in Figs. Gib) and^c), indicates very good agree- 
ment with that obtained via the imaging technique. This 
shape parametrization is not unique; an essentially indis- 
tinguishable source was also obtained for a fit performed 
with a Gaussian plus modified exponential 23] shape. 
The simpler Spheroid parametrization was chosen for the 
discussion below. 

The Spheroidal source function (Eq. [2Jl can be approx- 
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FIG. 2: (color online) (upper) Correlation functions, C(q), 
for 7r + 7r + and tt~it~ pairs and (lower) corresponding source 
functions, S(r) for (a, c) high kr most central collisions; (b, 
d) low kr peripheral collisions. Error bars indicate statistical 
errors with symbols as in Fig. 



imated by a short-range Gaussian source S sr (r); 



S sr (r) — A e 



/ {%Tt^R 2 T R a ) 



(3) 



for small r, and a long-range source Si r {r) for r » 
2R T R /^/R 2 - Rj given by 



S lr (r) ~ A ifee"' /4i? V (4^ (Jig - 



(4) 



Thus, the emission source shown in Figs, ^b) and^c) 
can be interpreted to reflect a short range Gaussian 
source of radius R sr = ^SR^R 2 / (2R 2 + R T ) and a long- 
range tail of extended space-time extent Ri r = Rq. The 
fraction of pairs associated with these sources X sr — 
A a 2 (3/(2a 2 + 1)) 3/2 and X lr = (A - A sr ), are obtained 
from Eq. |31 and the conservation of pairs respectively. 

Source functions were extracted via imaging and 
spheroid fits for several kx and centrality selections, in 
order to map the regions in kr and centrality where 
the long-range tail is prominent. Representative results 
are shown in Fig[3 for the indicated cuts. The experi- 
mental and restored correlation functions, compared in 
Figs[2Ja,b), indicate good agreement consistent within 
the statistics, as do the corresponding source functions 
shown in Figs|2fc,d). 

Figure |2| gives a more complete summary of the ex- 
tracted source parameters. The indicated systematic er- 
rors were obtained by varying the pair cuts in the analy- 
sis. The centrality and fcx dependence of the RMS radius 
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FIG. 3: (color online) (left) kr dependence of extracted 
spheroidal source parameters for pion source functions (a) Rt, 
(b) a, and (c) A for (filled circles) most central collisions; (open 
circles) perhipheral collisions, (right) centrality dependence 
of same parameters for (filled stars) low kr (0.2 — 0.4 GeV/c); 
(open stars) high kr (0.5 — 1.0 GeV/c). Filled bands indicate 
systematic errors. 



of the short-range source R sr , is similar (within 10%) to 
that obtained for i?i on g and i? s ido in an earlier analysis 
[l5| (cf. Figs. |3] (a) and (d)). The long-range source 
shows an effective radius Ri r — ax Rt, which is 3-4 
times Rt (cf. Figs. 0(b) and (e)) with the largest values 
for low kT and the most central collisions. The fraction 
of pairs exhibiting these characteristics is given by the A 
values shown in Figs. |3] (c, f); maximum prominence is 
shown for low fc^ pairs from central collisions. 

A central question is the origin of the long-range contri- 
bution to the emitting source. Instantaneous freeze-out 
of an isotropic source in the longitudinal co-moving sys- 
tem (LCMS) would give Ri r = R out = 7 x Rt in the 
PCMS. Thus, the values for a and 7, shown in Fig.^b), 
can be directly compared. At low kT, 7 is seen to be 
significantly less than a. Thus, a simple kinematic trans- 
formation from the LCMS to PCMS can not account for 
the observed source parameters at this kr- 

Could a composite particle emission source comprised 
of a central core and a halo of long-lived resonances ac- 
count for the observations? For such emissions, the pair- 
ing between pions from the core and secondary pions from 
the halo is expected to dominate the long-range emissions 
[24l | . If it is assumed that this halo includes only u> decay 
(ct ~ 24 fm), one may compare the tu yield with a simple 
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estimate of the fraction of pion pairs associated with the 
short- and long-range sources. Using the a and 7 values 
in Fig. Etb), A sr - 0.22 - 0.32, A ir - 0.23 - 0.13 
and Xi r /Xsr ~ 0-4 — 1-0 for the lowest kr- The pre- 
liminary data for lo/h~ =0.1 from Ref. |25j gives an 
estimate for Xi r /X sr ~ 2 x 0.1/\/A sr ~ 0.35 — .43 which 
lies at the lower end of the estimates obtained from the 
source parameters. Therefore, it is plausible that a max- 
imal kinematic boost combined with a halo of cos could 
account for Xi r /X sr . However, the steep centrality de- 
pendence of the radius of the long range source inferred 
from Figs. EJd,e) is not compatible with what would be 
expected for significant u> resonance contribution. Also, 
initial comparisons to the results from a recent dynamical 
model calculation |2y indicate that resonance emissions 
are not sufficient to fully account for the observed tails 
in the source functions. 

In summary, we have made the first extraction of the 
full ID emission source function for pions in Au+Au 
collisions at RHIC (^/snn = 200 GeV). This source 
function points to separate but prominent contributions 
from short-range emissions and a long-range tail of larger 
space-time extent than has been previously observed. 
This tail cannot be explained solely by simple kinematic 
considerations associated with a frame transformation 
from the LCMS to the PCMS. Further detailed 3D imag- 
ing measurements in conjunction with in depth model 
studies, are required to quantitatively pin down the ori- 
gin of this tail and to determine whether or not it is 
related to a phase transition. 



* Deceased 

' PHENIX Spokesperson zajc@nevis.columbia.edu 
[1] K. Adcox et al., Nucl. Phys. A757, 184 (2005). 
[2] E. V. Shuryak, Nucl. Phys. A750, 64 (2005). 
[3] S. Pratt, Phys. Rev. Lett. 53, 1219 (1984). 
[4] D. Teaney et al, Phys. Rev. Lett. 86, 4783 (2001). 
[5] T. Csorgo et al., Phys. Lett. B333, 494 (1994). 
[6] T. Csorgo, S. Hegyi, T. Novak, and W. A. Zajc (2005), 

nucl-th/0512060. 
[7] E. V. Shuryak, Phys. Lett. B44, 387 (1973). 
[8] U. Heinz et al., Ann. Rev. Nucl. Part. Sci. 49, 529 (1999). 
[9] T. Csorgo, Heavy Ion Phys. 15, 1 (2002). 
[10] M. A. Lisa et al., Annu. Rev. Nucl. Part. Sci. 55, 357 
(2005). 

[11] C. Adler et al., Phys. Rev. Lett. 87, 082301 (2001). 
[12] K. Adcox et al., Phys. Rev. Lett. 88, 192302 (2002). 
[13] J. Adams et al., Phys. Rev. Lett. 91, 262301 (2003). 
[14] Y. Sinyukov et al., Phys. Lett. B432, 248 (1998). 
[15] S. S. Adler et al., Phys. Rev. Lett. 93, 152302 (2004). 
[16] D. A. Brown et al., Phys. Lett. B398, 252 (1997). 
[17] D. A. Brown et al., Phys. Rev. C64, 014902 (2001). 
[18] S. Y. Panitkin et al., Phys. Rev. Lett. 87, 112304 (2001). 
[19] P. Chung et al, Phys. Rev. Lett. 91, 162301 (2003). 
[20] K. Adcox et al., Nucl. Instrum. Meth. A499, 469 (2003). 
[21] K. Adcox et al., Phys. Rev. C69, 024904 (2004). 
[22] S. E. Koonin, Phys. Lett. B70, 43 (1977). 
[23] P. Chung et al., Nucl. Phys. A749, 275 (2005). 
[24] S. Nickerson et al., Phys. Rev. C57, 3251 (1998). 
[25] W. Broniowski et al., Acta Phys. Hung. A22, 159 (2005). 
[26] A. Kisiel, W. Florkowski, and W. Broniowski (2006), 
nucl-th/0602039. 



